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Abstract: This note presents the Hilbert series technique to a wider audience in the context 
of constructing group-invariant Lagrangians. This technique provides a fast way to calculate 
the number of operators of a specified mass dimension for a given field content, and is a 
useful cross check on more well-known group theoretical methods. In addition, at least when 
restricted to invariants without derivatives, the Hilbert series technique supplies a robust way 
of counting invariants in scenarios which, due to the large number of fields involved or to 
high dimensional group representations, are intractable by traditional methods. We work out 
several practical examples. 
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1 Introduction 

In this post-Higgs-discovery era of fundamental physics, phenomenological models of physics 
beyond the Standard Model (BSM) are becoming increasingly baroque. The simplest models 
have been well-studied and more complex models offer numerous adjustable parameters that 
can be tuned to avoid ever more stringent experimental limits on new physics. One common 
way to augment the complexity of a model is to simply add particles transforming under 
higher-dimensional group representations. For some representative examples using sizable 
group representations in phenomenological settings, see [1-5] and references therein. In addi¬ 
tion to complex model building, another line of attack in the search for new physics takes a 
general bottom-up approach - parameterizing BSM effects through higher dimensional oper¬ 
ators formed from SM fields. These operators are suppressed by powers of some high energy 
scale, but they become important for detecting new physics in high precision measurements. 

Given a BSM model containing new particle multiplets in different group representations, 
one often wants the most general gauge-invariant Lagrangian containing all of the operators 
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allowed by the postulated symmetries. The formation of this Lagrangian allows a full explo¬ 
ration of the model’s experimental signatures. In the simplest cases, with small and familiar 
group representations such as SU{2) doublets, it is straight forward to assemble the La¬ 
grangian almost automatically. With larger group representations, the calculation becomes 
harder, and one often needs to take into account relations between group invariants (math¬ 
ematicians refer to these relations as syzygies) to obtain the correct number of independent 
terms in the Lagrangian. A similar level of calculational difficulty occurs in the task of form¬ 
ing higher dimensional operators from the Standard Model degrees of freedom. Even though 
the SM group representations are familiar, the high multiplicity of these representations con¬ 
tained in the operators quickly becomes challenging to deal with as the operator dimension 
increases. 

A fundamental reason for the prevalence of such group-theoretical calculations in phe¬ 
nomenological studies is the modern Wilsonian perspective of effective field theory. From 
this perspective, the gauge and global symmetries of the proposed field content completely 
determine the Lagrangian, and in principle one must include all possible invariant terms in 
every operator dimension as part of the resulting QFT. 

The mathematical tool known as the Hilbert series is perfectly suited for such computa¬ 
tions. The Hilbert series (or Molien or Poincare function) is a generating function encoding 
information about the number of independent group invariants that can be formed from some 
set of multiplets in different representations. As we explain in this paper, the Hilbert series 
provides an easy cross-check on calculations performed using more familiar group theoretical 
techniques, allowing one to ensure the correct number of independent terms in the Lagrangian. 
This technique is especially convenient when dealing with large representations or higher di¬ 
mensional operators, since the Hilbert series calculation is easily automated using computer 
algebra programs such as Mathematic A. 

The Hilbert series approach has been developed and used extensively in more formal 
theoretical settings. For example, it is often used in calculations involving the operator 
spectra of supersymmetric gauge theories [6-11], SUSY theories on D-branes [12-14], and 
moduli spaces of instantons or vortices [15-19]. Other references from a more mathematical 
point of view include [20, 21]. On the more phenomenological side, Hilbert series methods 
have been used to calculate the number of independent flavor invariants in the Standard Model 
and various extensions. The Hilbert series for leptonic flavor invariants in the Standard Model 
extended by the dimension-5 Weinberg operator was determined for two and three generations 
in [22]. This paper also obtained the Hilbert series of flavor invariants for the full type-I seesaw 
model in the case of two generations, and the series for quark sector flavor invariants for both 
two and three generations. This paper was followed by [23], which completed the calculation 
of the Hilbert series for the seesaw model with three generations, and also computed the series 
for quark flavor invariants with four quark generations. 

While the Hilbert series technique has been used in formal studies and in the more phe¬ 
nomenological setting of calculating flavor invariants, it is not yet well known and appreciated 
as a tool that can aid in the construction of general gauge-invariant Lagrangians. In this pa- 
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per, we aim to rectify this and to add the Hilbert series to the toolbox of a wider audience. 
Specifically, in Section 2 we introduce the basic Hilbert series concepts through a simple ex¬ 
ample, then present the general framework. Then, in Section 3, we work through a complete 
example to solidify the concepts. Sections 4 and 5 show some more complicated examples, 
and we discuss incorporating derivatives and equations of motion in Section 6. In Section 7, 
we conclude. Some mathematical background material is included in Appendix A. 


2 Hilbert series basics: [7(1) symmetry 


In this section we introduce the key mathematical ingredients of Hilbert series using a simple 
example. More formal introductory material to the Hilbert series can be found in the literature 
[7, 9, 10, 13, 21, 22]. 

For our example, consider a single complex scalar field charged under a 17(1) symmetry: 
0 —>■ —>■ The gauge invariant combinations are and there is exactly one 

possibility for each n. Writing the set of invariants as a series 

OO 

= ( 2 . 1 ) 

n=l 

where Cn is the number of different invariant possibilities for a given dimension, we have 

^ = 1 + (#*) + + (#*)3 + ... ( 2 . 2 ) 


Formally treating (</></>*) as numbers less than one^, 


this geometric series can be summed. 


= 


1 

1 — 


(2.3) 


Let us massage this further; the sum above can be replaced by an integral over 9, the variable 
that parameterizes the 17(1) transformation: 




2w 

j_ r _ <w _ 

2tt J {1 — (j) — (j)* e~^^) 

0 


(2.4) 


Substituting z = e*®, the d9 integral becomes a contour integral around \z\ = 1. 


i3 = 


1 


dz 


1 




(2.5) 


^To clarify that the objects we are manipulating, e.g. cf), are complex numbers rather than quantum fields, 
we will refer to (f) and similar objects as spurious. 
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The second piece of the denominator in Eq. (2.5) can be replaced by 


1 

(1 - (pz){l 


^ = exp I - log (1 - (j)z) - log (1 - y )| 



( 2 . 6 ) 


To get a better idea of what’s going on, let us expand the LHS of Eq. (2.6), again treating cf) 
and (f)* as small, complex numbers rather than quantum fields. To cubic order in both (j) and 


--—-^ = (1 + + (#*)^ + H- )+z (</> + </) (#*) + (t> H—) 

{l-(t)z){l - 

+ Z ^ ((/.2 + (/.2 (</,,/,*) + ...)+ ,/.3 + ^ + 1 (())*2 + <^*2 ^ ^ 

+ ^{(j)* + (t>* (#*) + (t>* + •••)■ (2-7) 

The series of U{1) invariants sits in the term with no 2 : factors, hence it is picked out when 
we multiply by l/z and perform the contour integral in dz (Eq. (2.5)). However, inspection 
of the expansion of Eq. (2.6) shows it contains all possible arrangements of 4> and 4>*; if we 
wanted to pick out the series of charge +1 combinations, we would simply have to multiply 
by ^ before taking the contour integral, the charge —2 can be accessed by multiplying by 2 , 
etc. 

In picking out a particular charge from Eq. (2.6), we are using two different mathematical 
facts. The first, already mentioned, is that Eq. (2.6) generates all possible combinations of cj) 
and (j)*, organized by charge. This exponential form of Eq. (2.6) is an example of a plethystic 
exponential [7, 9, 10, 13, 21], a generating function of all symmetric combinations of its 
argument (see A. 2). The second mathematical construction we employ is integration over the 
group volume, d9 = ^ for 17(1). Integrated over d9, terms containing any non-trivial power 
of 2 ; —>■ become integrals d9 e^^^ for some integer n and are therefore zero. Terms with no 
powers of 2 ; - the 17(1) invariants - remain and are the Hilbert series Sj. 

While the series of invariants in the 17(1) example above could be found without the aid 
of Eqs. (2.5) and (2.6), the power of this approach lies in its generality. Both the plethystic 
exponential and the integral over the group parameter 9 can be extended to sets of spurions 
transforming under arbitrary representations of arbitrary compact Lie groups. In addition to 
generating all possible combinations of spurions such as (j) and (p*, the plethystic exponential 
keeps track of relations among invariants, or syzygies (see Ref. [22] for some simple examples of 
relations among invariants). Our U{1) example is too simplistic to see these syzygies, however 
we will run into relations among invariants when we consider more complicated setups as in 
Sections 4 and 5. 

To create the plethystic exponential (PE) for a spurion A transforming with representa- 
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tion i? of a connected Lie group, we take 


PE[^, R] = exp (f; ^ ■ (2.8) 

r=l 

Here XRi^j) is the character of the representation R expanded as a monomial function of the 
j complex variables on the Cartan sub-algebra (equivalent to the group rank). For example, 
consider A to be in the fundamental representation of SU{2), a rank-1 group. The character 
for the fundamental representation is 

{z + l), (2.9) 

where z is a complex number with modulus one (often called a fugacity). The argument of 
the plethystic exponential is then 


E 

r=l 




z"' + 


+ =-log(l-^A)-log(l-H^). 


( 2 . 10 ) 


r=l 


This can be easily extended to more spurious and representations of different groups - all we 
need is a list of the characters for different representations of Lie groups. For a brief review 
of characters and further discussion of the PE, see Appendix A. 

Next, the invariants from the PE are picked out by the fact that the characters of compact 
Lie groups form an orthonormal set of basis functions on the Cartan sub-algebra variables 
[24] . As such, any function / of the sub-algebra variables can be expanded in terms of them: 

/(^i) = j= 6ij, (2.11) 


where the Cj are coefficients. The integration dfi in Eq. (2.11) is over the Haar measure - the 
volume of the group in question projected onto the Cartan sub-algebra variables (maximal 
torus). These volume elements can be found, for example, in Ref. [10] and are included in 
Appendix A.3 for convenience. While the PE contains all possible tensor products of spurious, 
the group integration projects out only the invariant combinations, resulting in the Hilbert 
series Sj. 

Looking back at our 17(1) example, we can rephrase the results in this more general 
language. The character for a 17(1) representation with charge Q is z^, so the PE in Eq. (2.6) 
is the sum of a representation with charge Q = +1 and representation Q = —1. The Haar 
measure for 17(1) is dzjz^ and the character orthogonality relation is the usual Fourier series 
orthogonality, as is best seen by setting z = e*®: 


27r 

J ^e^HQ-Q') =s{Q-Q'). 

0=0 


( 2 . 12 ) 
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Continuing, Eq. (2.5) can be understood as performing a Fourier expansion of the PE (Eq. (2.6)), 
multiplying by the trivial representation (i.e. 1) and integrating over the group. Since the 
characters are orthonormal, only the part of the PE expansion that lies in the trivial represen¬ 
tation ~ the 17(1) invariants - is projected out. Had we wanted to project out a different part 
of the PE, all we have to do is multiply by the conjugate character before doing the group 
integration. For example, to project out the charge -|-1 combinations, we need to multiply by 
Ijz = before integrating. 

Having shown how the PE and group integration can generate the Hilbert series of invari¬ 
ants for this simple 17(1) example, we now want to apply this method to the Standard Model 
(and its extensions). Specifically, we will take some subset of the fields Q,u^,d'^, L,e‘^, H, 
etc. of the SM (plus any extensions) as the spurions, dress them with the characters appro¬ 
priate to their groups and representations, then form the PE and do the contour integrations. 
The different SM spurions can transform differently, and some may transform under multiple 
groups. A separate group integration is carried out for each group under which the fields 
transform. If the contour integrals can be directly calculated, the resulting Hilbert series 
will be a rational function of the input spurions, ■ ■). This function 

S){Q, d‘^, L, ■ ■) can then be Taylor expanded as a multivariate power series in the 

spurions. The coefficient of any particular combination of spurions gives the number of in¬ 
dependent group-singlet operators that can be formed from the fields represented by those 
spurions. For example, a term 2Ql QL^L in the Hilbert series indicates there are two in¬ 
dependent singlets that can be formed from one Q field, one L field, and their hermitean 
conjugates. We emphasize that, while the Hilbert series gives the number of invariants, it 
does not give the particular index structure, so that must be worked out separately. Before 
we dive in, a few comments are in order: 


• Several of the spurions we want to use are fermionic, meaning they represent fermion 
fields. To properly count invariants including fermions, we need to extend the PE to 
handle antisymmetric spurions. This can be done by using the fermionic plethystic 
exponential (PEE) [21]. For a fermionic spurion A in representation R: 


OO 

PEF[A,R] =exp(^^ 

r=l 


{-1Y+^A^Xr{z^,) \ 

r ) ’ 


(2.13) 


where, as before, xnizj) is the character for representation 72 as a function of the Cartan 
sub-algebra variables. 

• As we are working with fermions (and, eventually, field-strength tensors) we must in¬ 
clude the Lorentz group representations for these spurions. At first sight, this seems 
problematic since the orthonormality of the group characters (the Peter-Weyl theorem) 
only holds for compact Lie groups. However, since our purpose is solely to count in¬ 
variants and does not involve any dynamics, we can work in Euclidean space, where 
the Lorentz group is compact: S'0(4) = SU{2)ji x SU{2)l- For simplicity we will take 


- 6 - 



all fermions to be left-handed objects, transforming as (0, under SU{2)l x 5?7(2)r; 
hermitian conjugate fermions therefore transform in the right-handed (^,0) represen¬ 
tation. When considering field strength tensors, we will work with the objects = 
X^i, ± iX^jy, which transform in the (1,0) and (0,1) representations. 

• While the PE or PEF generate all possible combinations of spurious, these constructs 
contain no information regarding the equations of motion - operator relations that go 
beyond symmetries. How to include the equations of motion, and derivatives in general, 
lies beyond the scope of this paper but will appear in a companion paper [25]. A few 
comments will be sketched out in the discussion in Section 6. 


3 Using the Hilbert series: a toy example 


It is instructive to work through a full example of the Hilbert series technique in order 
to demonstrate the general procedure. For this purpose, consider the Standard Model left- 
handed fermion doublet L, which is a weak isodoublet and a color singlet. In this example, we 
will ignore hypercharge, and we denote the number of generations as Xf. Since L transforms in 
the fundamental representation of both SU{2)l and SU{2)\y, the argument of the plethystic 
exponential is 

where x is the complex variable for SU{2)l and y is the variable for SU{2)\y- Recall that 
when carrying out the sum for the PE, we need to include the factor (—since L is 
fermionic. Explicitly, using the Haar measure for the SU{2) groups as given in Eq. (A.7), the 
Hilbert series for the spurion L is 


Sj{L) = 


(2vri)2 


J\x\ = l ^ J\y\ = l y 


exp 


"rE 


(- 1 ) 


r+l 


r=l 


U x’’ + 


2/^ + 


1 


X' / V y 

(3.2) 

Note that Xf is not a variable in the PE, but is rather a free parameter. Expanding the PE 
gives 


exp 




(-1) 


r+l 


r=l 


lM x^ + 


1++ 

xyj 


y ) 


Nf 


(^l + ^y\l + Lxyf^. 


(3.3) 


In general, it can be computationally challenging to do the contour integrals for the 
Haar integration directly, especially if a large number of fields is under consideration. In this 
situation, it is better to first expand the integrand in a Taylor series in the spurion L, and 
then integrate term-by-term up to the desired power [6]. The poles of the integrand will now 
all be at X = 0 and y = 0, so the residue calculation is much easier. This method will not 
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give the complete generating function ij, but rather the first few terms in its series expansion. 
This is sufficient for most purposes. However, if some all-order information about the series 
is needed, or the asymptotic form of the coefficients as the expansion variables go to infinity 
is desired,^ it will be necessary to directly do the integrals and obtain the full functional form 
of Sj. 

For the specific example of Eq. (3.2), we can illustrate this expansion method by leaving 
A^f unspecified and calculating the first few terms of the Hilbert series for L: 

Sj(L) = ^(Nf^-N{)L^ + ^(Nf^+5Nf^)L^+j^(5Nf^-3Nf^+17Nf^-69Nf^+50Nf^)L^+... (3.4) 

If A^f is specified, the Hilbert series can easily be calculated exactly, since the only poles in 
the integrand of Eq. (3.2) are at x = 0 and y = 0. Carrying out the contour integrations for 
the first few values of Nf gives the results displayed in Table 1. 


Af 

Hilbert series S){L) 

1 

1 

2 

l + L'^ + QL^ + L^ + 

3 

1 -b 3L2 + 21L^ -b 20L6 -b 21L8 -b 3L^° -b 

4 

1 -b 6L2 -b 56L^ -b 126L'5 -b 210L® -b 126^^° -b 56L^2 g^i4 ^16 


Table 1. The Hilbert series for the first few values of JVf, calculated from Eq. (3.2). Note that the 
resulting polynomials are palindromic, as noted in [22]. Also, for a given IVf there is a maximum 
dimension an operator can have. 


The entries in Table 1 can be understood using standard group theory. In this discussion, 
5C/(2)w indices will be explicitly displayed, and SU(2)l (Lorentz group) indices will not be 
displayed, but will be contracted within parentheses. For N{ = 1, there is no possible 
operator, since eai 3 {L'^Ll^) = 0. At the L‘^ level, the indices can be contracted in a single 
non-zero manner: ea 0 ^ 5 -y{L°^L^){L^L"^)- All other methods of index contraction for are 
either automatically zero, or can be written as this one by using the Schouten identities from 
the two SU{2) groups. At higher orders, for example L^, it is always possible to use the 
freedom given by the two sets of Schouten identities to rewrite the contractions in a form 
that is manifestly identically zero, which is why the Hilbert series terminates at order L^. 

For Af = 2, there is a single operator of order L^, namely eapiL^L^)- The fact that there 
are two distinct fields keeps this from being identically zero. Similarly, at order L^, the six 

^To obtain asymptotic information about generating functions, see for example the methods described in 
[ 26 ]. 









possible operators are 


)(L]L^), 

e„^e^,(L?Lf)(L7L^), )(^Ti2), ea^eps{L^Ll){LjLi). 

As a final illustration, for Af = 3, the three operators of order are ea/siL^L^), eap{L^L^), 
and ea 0 {L^L^)- 

The Hilbert series method can actually go further and give us the exact flavor content of 
the operators as in Eq. (3.5) directly, by doing something called “refining” the series. To do 
this, instead of putting a generic spurion L into the PE and multiplying by Af as in Eq. (3.2), 
we put in Af distinct spurious with different labels. Eor example, if Af = 3, use spurions Li, 
L 2 , and L 3 in the PE. This is also referred to as a “multigraded” Hilbert series. Doing this 
will in general increase the complexity of the residues that must be computed, so in some 
cases it may make the problem intractable. However, for the series of Eq. (3.2) it is feasible, 
at least for low values of Af. Calculating the multigraded series for Af = 2 and picking out 
the terms replicates the results of Eq. (3.5). Of course, refining the series does not give 
the exact index contraction structure (since there are multiple equivalent ways to contract 
the indices), but it does tell us that, for example, there are two independent operators with 
the flavor content L\L‘^. 

Because of the increase in computational power required by refining the Hilbert series, it 
may sometimes be advantageous to go in the other direction and “unrefine” the series. For 
example, if the Hilbert series involves both the lepton doublet L and the quark doublet Q, it 
can be useful to not only lump all of the flavors together, but also to lump together L and 
Q by setting both spurions to a common label t in the PE. This can significantly decrease 
the computational time required by a progam like Mathematic A, and may give all of the 
information that is necessary for a specific application of the Hilbert series method. 

The example worked out in this section shows that while the results obtained from the 
Hilbert series can indeed be replicated using other techniques, the calculations are often 
tedious and it is easy to make mistakes and miscount the operators. For example, showing 
that there are 126 independent operators of order for Af = 4 would be a formidable 
task by standard methods. By refining the series, the Hilbert series method provides an easy 
automatic way of obtaining the desired operators. 

4 Using Hilbert series for Standard Model effective field theory 

The Standard Model effective field theory (SMEFT) consists of the Standard Model La- 
grangian CgM plus operators with (mass) dimension greater than four that are invariant un¬ 
der the Standard Model gauge group SU{3)c ® SU{ 2 )w (g) U{1)y and contain only Standard 
Model degrees of freedom. At dimension 5, there is only one possible operator, the Weinberg 
neutrino-mass operator [27]. Continuing on to dimension 6 , a classification of the available 
operators was done in [28], and the reduction to a minimal set of dimension -6 operators was 
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carried out in [29], resulting in a set of 63 operators. More recently, the construction of the 
minimal set of dimension-7 operators was also completed, giving 20 independent operators 
[30]. Any BSM physics can be matched onto this effective field theory by integrating out the 
postulated new heavy particles. In the past few years, much work has been done towards 
understanding the structure and use of the SMEFT [31-41]. 

Calculating the full set of operators for a given dimension with the Hilbert series requires 
including covariant derivatives, field strengths, and the equations of motion. Some thoughts 
on these ingredients are included in Section 6, and we leave an in-depth discussion to future 
work [25]. For simplicity, in this section we focus on operators without these complicating 
factors. This still provides a good example of the practical applications of the technique. 

4.1 Dimension-6 baryon number violating operators 

The Hilbert series technique proves to be useful for finding independent SMEFT opera¬ 
tors, even when restricted to a subset of the field content. For example, consider the set 
of dimension-6 baryon-number-violating operators. Ref. [29] presented five such operators. 
However, it has been noted that only four of these five operators are independent if the flavor 
structure is taken into account [37, 42]. It would be nice to see such operator relations without 
doing a detailed calculation, and indeed using the Hilbert series can bring the dependence 
betwen different operator structures to light. 

Consider the class of baryon-number violating dimension-6 operators with field content 
QQQL, where Q is the left-handed quark doublet and L is the lepton doublet. Taking Af = 3, 
the argument of the PE for these fields is 

^+ f + i) (x + (y + (4.1) 

where x is the variable for SU{2)w, y is the variable for SU{2 )l, u is the variable for t/(l)y, 
and { 21 , 22 } are the variables for SU{3)c- Calculating the multigraded Hilbert series so that 
we can see the form of the operators gives 

^{L, Q) = 1 + 57 LQ^ + 4818 L‘^Q^ + 162774 + ..., (4.2) 

so we expect 57 independent operators of the form QQQL when flavor structure is included. 
Now consider the specihc operator structure 

■ (4.3) 

This is symmetric in the flavor indices {p,r}, so without even considering additional sym¬ 
metries implied by SU{2)l Fierz identities, the maximum number of independent flavor per¬ 
mutations is 3 X 3 X 6 = 54 < 57 (since a 3 x 3 symmetric matrix has only 6 independent 
entries). Therefore the operator structure in Eq. (4.3) does not capture all of the 57 inde¬ 
pendent operators that we know exist at this order. Perhaps changing the SU{2)w index 
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contraction structure will enable us to capture all of the operators in a single expression. To 
this end, consider the following structure with the 5?7(2)^ contractions “offset” from the 
SU{2)l contractions: 

eaf^-yeike,i . (4.4) 

None of the symmetries are immediately apparent in this form, so naively one might think 
that we can now just multiply out the flavor possibilities and get 3^ = 81 operators with 
different flavor structures. But this neglects SU{2)l Fierz identities (see for example [43]), 
the relevant one of which gives 






(4.5) 


Rearranging the first term on the right-hand side using the SU (2)w Schouten identity tijemn = 
^im^jn — ^in^jrm and relabelling the second term on the right at the cost of a sign gives 






2 (g;“gf) (g^4) - (grgf) (g^^^) 


(4.6) 

This identity tells us that among the 81 different flavor permutations of the structure in 
Eq. (4.4), there are 3^ = 27 relations among operators differing only in Q flavor permutations. 
When all of the quark flavors are identical (p = s = r), the relation is trivial. Since (p = s = r) 
occurs 3 times for iVf = 3 as in the Standard Model, there are 27 — 3 = 24 linear relations 
among the 81 flavor permutations, leaving 81 — 24 = 57 independent permutations! Since 
the Hilbert series in Eq. (4.2) revealed that there are only 57 independent operators, we are 
done; the structure in Eq. (4.4) contains all of the possibilities, and we do not need to write 
down a different fermion current structure or different SU{2)\y index contractions. 

The Hilbert series technique is by no means restricted to invariants of gauge symmetries; it 
can deal just as easily with global symmetries. As one example, we could have included baryon 
number for all fields, adding a character and group integration for that global U{1). Had we 
included baryon number in the above example, we would find no invariants - by construction, 
since our example concerned baryon number violation. However, had we considered a wider set 
of spurions, both baryon number violating and respecting terms in the Hilbert series would be 
generated. In this case, integration over baryon number could be used to project out different 
subsets, e.g. only baryon number respecting operators, operators violating baryon number by 
one unit, by two units, etc. Of course, the same technique could be used for finding operators 
violating lepton number or B — L. 


4.2 The dimension-7 operator LLLeH 

As another example of applying the Hilbert series to the SMEFT, consider the class of 
dimension-7 operators with field content LLLeH [30]. Here e is a left-handed field which 
can also be written as and H is the Higgs doublet. In the argument of the plethystic 
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exponential we then have (again with Nf = 3): 

3L + 3e tt , (4.7) 

where x is the variable for SU{2)w, y is the variable for SU{2)i, and u is the variable for 
C/(l)y. Calculating the unrefined Hilbert series gives 

f)(t) = 1 + 3 + 57 + 171 + 6 + 144t^ + 1053 + ... (4.8) 

Order contains the dimension-7 operators, and calculating the multi-graded Hilbert series 
shows that this order contains only the dimension-7 operators of the form L^eH, so we know 
that there are 57 operators with various flavor structures. 

One such possible structure is 

eijemn {epLi) (L^,L™) (4.9) 

As in the previous section, the symmetries are not manifest, so we need to use a Fierz identity: 


€7 4 6- 


ij ^mn 


(epLi) {LiLf) = -eijemn {epL{) {LfLl) - eijemn {epLf) [V,L{) (4.10) 


In the first term on the right-hand side of Eq. (4.10), we can switch the {i,j} labels at the 
cost of a sign (and use the fact that {Z 1 Z 2 ) = (^ 2 - 21 ) for anticommuting 2-component spinors). 
The second term on the right-hand side can be rewritten in a similar way after applying the 
SU{2)]v Schouten identity. The final result is 


CijC 


tj ^mn \ 


iepL\){Li.LT)H^ = [{epL\){LiLT) + {epLl){LiL^) - {epL\){L{L^)] . 

(4.11) 

In a similar way as in the previous section, this gives 27 relations among the 81 flavor permu¬ 
tations, this time among operators differing in L flavor permutations. Again only 24 relations 
are non-trivial, leaving 57 independent operators. So the structure in Eq. (4.9) encapsulates 
all of the independent operators of the form LLLeH. 

Now suppose that instead of starting with the structure in Eq. (4.9), the first guess had 
been instead 

eijemn{epL^){KLi)H^. (4.12) 


Without even considering Eierz identities, we can see that this structure is antisymmetric in 
{r, s}, so it contains a maximum of 3 x 3 x 3 = 27 independent flavor permutations, since 
a 3 X 3 antisymmetric matrix has only 3 independent entries. Thus we know that another 
structure is needed in order to get the full 57 independent flavor permutations. In this way 
the Hilbert series allows a check on the generality of a specific SU{2)l and SU{2)w index 
contraction structure. 
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5 Hilbert Series for BSM 


In the phenomenological study of extensions to the Higgs sector, it is often necessary to write 
down the most general form of the Higgs potential including various multiplets of SU{2)\y 
with potentially different hypercharges. Usually only renormalizable terms are included in 
the potential, but in some cases higher-dimensional terms are also necessary. In either case, 
this exercise is a perfect candidate for using the Hilbert series to check that the Lagrangian 
includes a complete set of operators up to a given dimension. 

As an explicit example, consider extending the SM Higgs sector by adding a scalar mul- 
tiplet X which is a quadruplet under SU{2)\y and has hypercharge —1/2, as was done in [1]. 
The dimension-2 terms for the Higgs sector potential are then trivial to construct, and there 
are no possible dimension-3 terms, but the task becomes more complicated at dimension 4. 
Using the character function for the 517(2) quadruplet as given in Appendix A (Table 2), the 
argument of the PE for the SM Higgs <I> and the new field x is 


X{z-^ + z + - + 






u 


1/2 _ 


+$ u+- 


Generating the first few terms of the unrefined Hilbert series gives 


,-V2_ 

(5.1) 


Sj{t) = l + 2t‘^ + llt^+ 31t^ + 94t^+ 222 t^^ + ..., (5.2) 


so we see that there are 11 independent operators at the dimension-4 level. ^ 

As a side note, the field and group representation content for this example is simple 
enough that the closed form for the unrefined Hilbert series can be calculated. The result is 


Sj{t) 


l + -h 17t® -h 13t^° -h 17t^2 9^14 4^16 ^20 

(t4+t2 + l)2(i2 + l)4(i2_l)8 


(5.3) 


Since the complete Hilbert series can be calculated, and since the result is a rational function 
of t, it is always possible to calculate a closed form for the coefficients of the series (see 
for example Chapter 7 of [44]) . One way to do this is to use the Mathematic A function 
SeriesCoeff icient, putting in a general integer n as the argument for the coefficient order. 
The resulting expression for the coefficients of Eq. (5.3) is complicated and not that useful, 
since calculating the Taylor series to the necessary order is simple in this case. However, 
there might be situations where calculating a closed form for the coefficients could be useful 
for figuring out general information about the form of the series such as the asymptotic form 
of the coefficients. 

Returning to the calculation of the (renormalizable) Higgs potential, by using Mathe- 
MATICA or a similar program, we can generate the multi-graded Hilbert series and extract 

®As all operators are scalars, counting dimensions is the same as counting the number of spurious. This 
can be generalized by weighting each spurion by its mass dimension, i.e H eH, Q —l then collecting 

like terms. 
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the relevant terms, which gives operators with the following field content: 

2xx^xx^ ^'^XX^X^ d>^$^X^X^ ^XXX^ 2d>d>'^xx^ 

<h<i>xX) <i><h<i)l^<i>l^. 

Counting up the terms we do indeed get eleven. Note that for the field content XX^XX^ the 
Hilbert series tells us that there are two independent terms, which is indeed the conclusion 
reached in the paragraph following Eq. (13) of Ref. [1]. The same conclusion holds for the 
field content d>t<f>x^X- The Hilbert series provides these results without the need for any 
detailed calculations involving group invariants. These calculations can of course be done 
without much trouble for these dimension-4 operators, but it is nice to have this check. Also 
note that the Hilbert series result (5.4) gives an operator and its hermitian conjugate that 
are not included in the result for the most general renormalizable Higgs potential given in 
Eq. (12) of Ref. [1], namely the operator with field content XXX^^- Working out the correct 
SU{2)w index contractions necessary (using the same index notation as Ref. [1]) gives the 
operator structure 

\.ijkXk'mnXijm^^'^nn' T h.C. (5'5) 

Using the Hilbert series has allowed the quick identification of an operator that was missed 
in the original potential; a clear demonstration of the practical efficacy of this technique. 

The 31 independent dimension-6 operators for this example can also be easily generated. 
The resulting operators are 

3(x^x)^ 2 (4>t)2x(x^)^ 2$x^(x^)^ 

4<h4.tx2(x^)2, 3 cl>(cl>t)2;^(xt)2, ci)(<ht)3(xt)2, 2<h2xV, 3<h2chtxV, (5.6) 

6 Thoughts on derivatives 

The examples we have worked through so far have been limited to collections of spurions 
without derivatives. Derivatives are necessary if we wish to apply the Hilbert series to a 
wider set of problems, such as the full set of SHEET operators of a given dimension. An 
immediate complication when including derivatives are the equations of motion (EOM) - 
relations among operators that are not governed by symmetries or invariances. 

At first glance, it seems like the Hilbert series is ill-equipped to handle EOM. However, 
the Hilbert series is not a Lagrangian and is merely a tool to count invariants formed from 
whatever spurions are put into the PE/PEF. For the purposes of counting invariants, the 
role of the EOM is to remove spurions namely derivatives on fields - by swapping them 
for different combinations of spurions with no derivatives, which is something that can be 
handled by the Hilbert series methodology. Consider the EOM for the left-handed Standard 
Model quark doublet Q: 


ilpQ = ylu‘^UH* + yld<^^H, 


(6.1) 


which allows one to remove lJ)Q from the set of spurions. Derivatives of fermions can still 
appear in invariants, first showing up at dimension 7 [30]. For example, consider the operator 

OLLud = eij ( 6 . 2 ) 

The difference between the derivatives in the two equations above lies in the Lorentz group 
representation. Acting on a left-handed fermion with a derivative, we get 



The EOM in Eq. (6.1) involves the (^,0) representation only. Therefore we can incorporate 
the EOM for a fermion field il) by including an additional spurion for the (^, 1) part of 
but omitting the (^, 0 ) spurion. 

If we follow this logic, as we add further derivatives etc. we should add a new 

spurion to the PE/PEE for every new representation of the Lorentz group that is formed. 
Eor derivatives of scalars such as the Higgs field H, D^H is new and should be added to the 
PE. Eurther derivatives, such as will contain a piece with (0,0) Lorentz structure; 

this piece should not be included as a separate spurion as it is eliminated by the Higgs EOM, 
□H = m?H -!-•••. Similarly, for derivatives of the field strength tensors the ( 5 , 5 ) 

representation should be omitted. 

Even after removing the derivative spurions, the Hilbert series still lacks information 
about integration by parts, so invariants involving derivatives need to be checked manually 
for redundancy. This check will get tedious if we Taylor expand the Hilbert series to arbitrarily 
high mass dimension operators, but it should be reasonably manageable for the dimension -8 
operators. Eurther investigation along these lines will appear in Ref. [25]. 

7 Conclusion 

The Hilbert series is a mathematical method providing a generating function for group in¬ 
variants. Since this method is not well known in the phenomenological community, we have 
introduced the Hilbert series through simple examples and provided several practical illus¬ 
trations of its use. Through the presented cases, we see that the Hilbert series technique 
proves very useful for common applications. Any calculation where a BSM or effective field 
theory Lagrangian contains fields transforming under uncommon group representations or 
where operators of higher order in mass dimension are needed, the Hilbert series is an easy 
way to get the right number of invariants. This is important in particular when looking for 
small deviations from known physics by comparing to precision measurements as redundant 
or missing terms in Lagrangians could invalidate analyses. 

With the inclusion of derivatives and the equations of motion, and applying sufficient 
skill towards reducing the computational load of calculating the necessary residues, it should 
be possible to compute a complete Hilbert series for the Standard Model effective field theory. 
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This generating function would contain at each dimension a coefficient specifying the total 
number of independent operators of that dimension. If the task was extended to calculate 
the multi-graded Hilbert series, one would then be able to directly see the different forms 
of operators at each order simply by expanding out the Hilbert series. We intend to work 
towards this goal in an upcoming study [25]. 
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A Mathematical background 
A.l Characters 

The character of an irreducible group representation is the trace of the matrix giving the 
representation. Any matrix representing a group element in a given representation will have 
the same character. 

In general, the characters of irreducible Lie group representations are obtained from the 
Gartan matrices, using the Freudenthal recursion formula to find the correct multiplicities of 
weights (see [21] and references therein). In practice, this calculation can be trivially done 
using the Mathematica package LieART by utilizing the WeightSystem command [45]. The 
results for some common representations are shown in Table 2. The characters can also be 
found using the character generating functions outlined in [ 21 ]. 


Representation 

Character function 

SU{2) fundamental 


SU{2) quadruplet 


SU{3) fundamental 

^i + t + i 

SU{3) anti-fundamental 

^ + ^2 

SU{3) adjoint 

"i"2 + | + S + 2 + | + | + ^ 

SU{2) adjoint 

z^ + l + j. 

1/(1) with charge Q 



Table 2. Character functions for several common group representations. Note that setting all of the 
variables to 1 gives the dimension of the representation. 
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Many common group calculations can be done using characters. For example, consider 
taking the tensor product of two SU{2) fundamentals. Using the characters from Table 2, 
this looks like 

which is just the familiar triplet-singlet spin decomposition ^ (8) ^ = 0 © 1. 


A.2 The plethystic exponential 

For a multivariable function /(U,..., tn) satisfying the property of going to zero at the origin, 
the plethystic exponential (PE) [7, 9, 10, 13, 21] is 


PE[/(ti,...,tn)] =exp 



(A.2) 


The plethystic exponential generates all symmetric combinations of the variables of the func¬ 
tion f{ti, ..., tn)- For example, the plethystic exponential of f{A, B) = A + B is 

+ ^] = (1 _ _ B) =1 + A + B + A^ + AB + B^ + ... (A. 

For fermionic variables we are interested in antisymmetric combinations instead of sym¬ 
metric combinations, since fermions obey Fermi-Dirac statistics. The fermionic plethystic 
exponential is defined as [21]: 



PEF[/(ti,.. .,tn)] = exp 


' OO , ^ 


(A.4) 


\r=l 


Eor example, 


PEE[A + B] = l + A + B + AB, 


(A.5) 


where the first three terms are trivial, and the last term is indeed antisymmetric under A ^ B, 
since fermions anticommute. In the body of this paper, the fermionic plethystic exponential 
is not always explicitly differentiated from the ordinary plethystic exponential, but the factor 
of (—1)^"*“^ is always included with fermionic variables (spurious). 


A.3 Group integration with the Haar measure 

It is possible to integrate over the manifold of a Lie group by using an invariant measure 
known as the Haar measure. This group integration projects out invariant quantities from 
the combinatorial expansion of characters provided by the plethystic exponential. The Haar 
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measures that are used in this paper are as follows (taken from [10]): 



(A.6) 

(A.7) 



d/^5C/(3) 


1 

(27ri)2 


dzi 

|2l| = l 



dZ2 

Z2 


(1 - Z1Z2) 



(A.8) 


Further examples of Haar measures for various Lie groups can be found in [10], where a 
general formula is also presented. 
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